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The' compulov proRmui HTGAP deacrlbod in this user's manual was pre- 
pared under NASA Contract NAS5-2A14^», "Axially Grooved Heat Pipe". The 
contract was administered by Goddard Space Flif^ht Center, Greenhelt, 
Maryland, with Mr. Stanford 01 .endorf serving as technical monitor, 

Dr. V. Kamotani was re.spon.sihle for the development of the analytical 
model and computer pro^^ram. 



1 . INTKOIHIOTION 


UTOfU' is ji ov>minitin' in‘Oj\r,vm that nvnuorloally prinUots tho sttvuly 
stato tomjv'raturo ilist. rllnit Ion InaiiU' an nxlally Rrooveil t\oat plpo wall 
for a nlvon j^oomotry aiul workln}; flvriii umlor various heat input and 

output modes. The pi\\ktram computes both evaporat and condenser film 
coel f icients. The program is able to handle botfy axisywmetrlc and non- 
axlsymmeiric heat transfer cases, Non^axisymmetric heal (ransl’er restilts 
either ttoro non»nnUorm input at the evapor.'itor or non-uniform heat removal 
Irom the cvuidenser, ot front both. The presence of a liquid pool in the 
condenser re\tiv>n under one"j:t condition also c.auses non“,txls>i;m«etvle 
heat trauster, and its effect on the pipe wall temperature distribution 
is included in the present pro^r.'\m. 

The hydrodynamic aspect of an axially jtrooved heat pipe is studied 
in fite droove Analysis Proy,ram (UAT) . The present thermal analysis pro- 
jiram asstunes tliat the dAP prop.ram for othor similar prop,v.amsl is rv\n 
first so that tl»e heat transport limit and optimum fluid charjye of tiu' 
heat pipe are known a priori. The perfonnance of an uuder-chavped heat 
pipe is not eonsldered in the present program. 

d. .AN.M.YS1S 

» . V P rob 1 em d t ,a t emeu t 

Figure 1 shows a typical axially ,tu‘ooved heat pipe. Eaci\ p,roove is 
desijtnod to operate as an Isol.nred liquid p.assapie. I'nder .steady st.ate 



conditions the rates of condensation and evaporation must be equal for each 
groove. For uniform heating at the evaporator and uniform cooling at the 
condenser (axisymmetric case) all the grooves are subject to identical 
boundary conditions so that the performance of the heat pipe can be deduced 
from the behavior of any one groove (assuming no excess liquid) . For 
non-uniform heating and cooling, however, the rates of condensation and 
evaporation of a groove may be different from those of .ther grooves so 
that it is necessary to study the behaviors of all the grooves in order to 
predict the performance of the heat pipe. 

The pipe wall temperature distribution is governed by a steady state 
heat conduction equation together with t.he exterior boundary condition 
(i.e., the way heat is applied to the evaporator and removed from the 
condenser) and the interior boundary condition (i.e., the rates of 
evaporation and condensation). Using cylindrical coordinates (Fig. 1) 
the governing equation is written as 

j. ^ -I- J- ^ _i. £I == o 

cM'^ >-<57' <PZ^ (1) 

The boundary condition describing the thermal interaction between 

the groove and the liquid in it, and between the liquid and the vapor is 

very complicated because of the complex geometry of the groove. To avoid 

an unnecessary complicated solution in Che present analysis, the inner 

surface thermal behavior is modelled as that shown in Fig. 2. An 

equivalent heat transfer coefficient ( heg^e or ) is used to describe 

C 

the thermal interation between a hypothetical inner surface (radius ) 
and the vapor. This approach was initially used by Schneider and 
Yovanovich [1]. The computations of andjjeg.c ate explained in the 
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following section. The region covered by a liquid pool is assumed to be 
thermally insulated. 

The exterior boundary condition is determined by a specified (by the 
user) heater-cooler combination. Three different boundary conditions can 
be specified in the regions covered by the heater and cooler as explained 
in section 2.3. The rest of the pipe wall surface is assumed to be 
thermally insulated, or the user can specify parasitic heat loss or gain. 
The heat conduction through the pipe end caps is neglected in the present 
analysis . 

Equation (1) together with the above mentioned boundary conditions is 
solved numerically using the successive-over-relaxation method. 

2.2 Problem Formulation 

To solve Eq. (1) numerically the heat pipe wall is divided into a 
number nodes as shown in Fig. 3. The cylindrical wall is divided into L 
nodes radially (r-direction) , M nodes (equal to the number of grooves) 
circumferentially (Y^ -direction) , and N nodes axially (z-direction) . An 
energy balance applied to every node yields a system of finite difference 
equations representing the heat conduction equation (1) . The following 
finite difference formulation follows Schneider and Yovanovich [1], 

For the (i,j,k) th node the finite difference equation is written as 


= c, 

-h 


Cs TiJ-ly A 

~h 

Cif. 


-f 

C 6 Td, / 


( 2 ) 
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where 


c 


U2. - 


cfY^ 


n r) 


cfy^ 


h t- cfYV2. 
Yi - eft /z 


r 

L5,L 


( d~rV2 ) — (h - <f t~/2 ) 
z 


(3) 


The heat transfer rate at the hypothetical inner surface (Fig. 2) of 
the pipe wall is determined by the equivalent heat transfer coefficient 
Ileg represents the total thermal conductance between the pipe inner sur- 
face (radius Ri) and the vapor. Then, the interior boundary condition is 
incorporated into the above formulation by putting 

/le£ t 

C, = — ^ 


^,v’ 2- 


heg 



for those nodes contacting the vapor. In the present analysis the vapor 
region is modeled as one node with a uniform temperature T^. 

If the pipe inner surface temperature is lower than T^, condensation 
takes place. In Ref. 2 Kamotani suggested a numerical method to compute 
the condenser film coefficient of a grooved heat pipe. However, the method 
is too complex to be included in the present program. After some cal- 
culations using the method under various conditions, it was found that 
the following expression for the equivalent heat transfer coefficient 
expresses well the computed as well as experimental values. 
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ITLRi • 02 4 - M. ^ 

Aw S 


( 5 ) 


If the pipe inner surface temperature is larger than T^, evaporation 
takes place. After studying several experimental and analytical investi- 
gations on the evaporator film coefficient of a grooved heat pipe 
(Ref. 3-7), the following expression for the equivalent heat transfer 
coefficient is suggested. 



N /f/ / 

.07 + 4 -^ 

J^rV S 


( 6 ) 


The boundary condition at the exterior surface of the heat pipe is 
discussed in the following section. 


2 . 3 Modes of Heat Input and Removal 

The exterior boundary condition is determined by the way heat is 
applied to the evaporator and removed from the condenser. Three typical 
boundary conditions are considered in the present analysis. They are shown 
schematically in Fig. 4, and discussed below. 


Type 1 Boundary Condition 

Total heat flux Q and vapor temperature T^ are given. This boundary 
condition is incorporated into Eq. (2) by adding a source (or sink) term 

s = ± f 

to the right hand side of Eq. (2) for those nodes contacting the heater 
or cooler, q is heat flux per unit area (Q divided by heating or cooling 
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area). + sign is for heat input, and - sign for heat removal. 

Type 2 Boundary Condition 

Total heat flux Q at the evaporator and condenser surface temperature 
T^ ^ are given. For those nodes contacting the heater the treatment is 
the same as above. The surface temperatures of those nodes contacting 


the cooler are set equal to T . The vapor temperature is computed 

w , c 


from the relation 


U ( Ts-c - Tv ) 


O 


l/n-rtjM. /OA^'jlcc. 


where Tsl. is the inner surface temperature. 


Tyep 3 Boundary Condition 

Both evaporator and condenser surface temperatures T^ ^ and T^ ^ are 

given. The specification of this boundary condition and the computation 

of T are the same as above. 

V 

If there is parasitic heat loss or gain from the heat pipe to the 
environment, it is accounted for by putting 








for the nodes on the pipe outside surface, is the heat transfer 

coefficient between the pipe and the environment, and T . is the ambient 

amb 


temperature. 



2.4 Puddle Effect Treatment 


OF POOK r'S'-'^'-v 

* ^ 


At the user's option the program investigates the effect of a liquid 
pool on the pipe wall temperature distribution in one-g and under an over- 
charged condition. The excess liquid forms a liquid pool at the bottom 
of the heat pipe, which influences the pipe performance (puddle effect). 

The hydrodynamic aspect of the puddle flow is studied in GAP. The present 
program considers mainly the thermal aspect. 

The free surface of the liquid pool is assumed to be flat in the 
present analysis. The assumption is not valid when the heat pipe is 
operating at a small tilt angle under a high heat load, or when the ratio 
f ^ A v/V (Bond number) is on the order of one or less. Therefore to avoid 
a complex combined hydrodynamic and thermal analysis, the puddle effect 
at zero tilt angle is not considered in the present program. 

Depending on a given tilt angle and amount of excess fluid, the puddle 
shape is classified into one of four shapes depicted in Fig. 5. To 
simplify the numerical analysis the following approximate relation between 
the puddle volume and the puddle depth at the condenser end is used, 

( :iRy) 

The pipe inner surface submerged in the pool is assumed to be thermally 
insulated. Numerically is set equal to zero for the submerged nodes on 
the pipe inner surface. At a given z location the submerged region is 
specified by the angle Q (Fig. 5). The angle 0 is related to the puddle 
depth H by the equation 

0 = ^ ^ J 
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2.5 Method of Solution 

A set of finite difference equations (2) together with the aforementioned 
boundary conditions are solved by successive-over-relation to determine 
the pipe wall temperature distribution. Recognizing the fact that the 
heat transfer in the pipe wall is predominantly in the radial direction, the 
temperature distribution in the r-direction is determined first to reduce 
the convergence time (Schneider and Yovanovich [1]), that is, at a given 
(Y^,z) l:<cation^ L equations in the r-direction are solved simultaneously. 

Once the temperature distribution in the pipe wall is determined, the 
program computes the evaporator and condenser over-all film coefficients 
which are defined as 


L ^ a 

2K R^Lc 

Notice that the coefficients are based on the vapor core diameter (2R^) 
and the difference between the vapor temperature and the average evaporator 
or condenser surface temperature. 

The program also computes the heat load of each groove. In case of 
non- ax i symmetric heat transfer each groove transports different amounts of 
heat so that all the grooves do not dry-out at the same time. The program 
checks whether the heat load of any of the grooves exceeds the heat 
transport limit of a single groove. To compute the single groove trans- 
port limit the user must specify the heat transport limit of the heat pipe 
in zero-g (or under zero tili^ without the puddle effect) and under 
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axisymmetric heat transfer conditions, which can be determined either from 
GAP or from similar hydrodynamic analyses. If the transport limit is 
not known, it can be set equal to an arbitrary large value, but in that 
case it is noted that the program does not predict partial dry-out con- 
dition. In zero-g or for horizontal operations without the puddle effect 
in one-g, the single groove transport limit is simply equal to the heat 
pipe transport limit divided by the total number of grooves. Under tilt 
conditions in one-g without the puddle effect the single groove transport 
limit is calculated as 




f f j ^ 


/V' 


A's'7 


(8) 


In the above relation the static height of the pipe must be speci- 

fied by the user. .)(s-T can be determined from GAP or from the relation 

X = 

When the puddle effect is included, the single groove transport limit is 
calculated as (Ref. 8) 


9 _ Sld!li±Jr^ f , _ ) 

^ XsT J 


( 9 ) 


where 


Le L(\ ^ ( Lp - LcJ L f S. L c 

^ ^// ~ ^ ^ Lp— Lc/^ ‘-i 

V 2’(^r ^ ^ c ) ^ Lp >■ L Ci -h c. 



If the heat load of any of the grooves exceeds the transport limit, 
the program writes-out a statement to inform the user about the partial 
dry-out condition. 

It is noted that the above analysis is based on the assumption that 
there is a groove-to-groove liquid communication in the condenser section. 
According to Kamotani [9] such condition most likely prevails in one-g 
as well as in zero-g applications. Therefore the present program does 
not consider the situations where there is no liquid communication 
between the grooves. 

3. PROGRAM DESCRIPTION 

3.1 General 

HTGAP's flow chart is given in Appendix A, and the program listing in 
Appendix B. The program is written in Fortran IV, and is designed to 
operate on the IBM 360 computer system. The field length of the program 
requires 60 K actual works, and the complications of the program requires 
typically .5 CPU minutes and .5 10 minutes. The CPU time varies depending 
mainly on the desired accuracy in the numerical iteration. The variables 
and constants used in HTGAP are listed at the beginning of the program. 
HTGAP uses metric units. 

HTGAP consists of a main program and a subprogram. The main program 
reads the input data, conducts numerical iterations, and outputs the 
final results. If required, the subprogram (PUDDLE) is called, which 
determines the region covered by the puddle. 



3.2 Input Descilption 

The HTGAP user must correctly specify the properties of the working 
fluid and the pipe material, the pipe and groove geometry, the selection 
of the type of heat input and removal mode^and the operational mode of the 
program. The set-up of the input data deck is shown in Table 1. In 
case of multiple runs these input cards specified in Table 1 (also in the 
program) must be repeated for eich additional run. The rest of the para- 
meters remain fixed. 

The physical properties of the working fluid and the pipe material 

should be evaluated at the vapor temperature (T^) for type 1 boundary 

condition, and at T for types 2 and 3 (T is not known a priori in 

’ w,c V 

these cases). In the latter cases a better accuracy is obtained by running 
the program twice, evaluating the physical properties at the calculated 
vapor temperature from the first run. 

Care'S must be taken in specifying the heating and cooling zones . 

Four numbers are used to specify the circumferential extent of the heat- 
ing (or cooling) zone. They are angles (in degrees) measured counter- 
clockwise from the pipe bottom. Some examples are shown in Fig. 6. If 
only two numbers are needed (examples 2 and 3) , the rest of the numbers 
should be set equal ta^an aribt^'ary number Ig j wgto' r than 500. 

3.3 Output Description 

The program output from the illustrative example (explained in the 
following section) is shown in Appendix C. On the first output page the 
program outputs the input data consisting of the names of the working fluid 



and the pipe material, the pipe and groove geometry and the heating and 
cooling mode. Those parameters are the same for all the cases in case pf 
multiple runs. The following two pages describe the heat pipe operating 
conditions, the physical properties and the computational results for one 
run. The computational results include the heat pipe surface temperature 
distribution, the heat load of each groove and the pipe performance 
characteristics (average evaporator and condenser temperature drops, film 
coefficients, partial dry-out conditions, etc). 

3 . 4 Sample Case 

To illustrate HTGAP, the following case is studied. An axially 
grooved aluminum heat pipe filled with ammonia is operating at 250 K in 
one-g. 


heat pipe geometry 

evaporator length 
adiabatic length 
condenser length 
outer diameter 
inner diameter 
vapor core radius 

groove geometry 

number of grooves 
groove depth 
average land width 


Le = .3048 m 
La = ,4572 m 
Lc = ,1524 m 
2Rout = .159E-1 m 
2Rin = .1088 E-1 m 
R = .45E-2 m 

V 


N = 27 

D = .108E-2 m 
S = .37E-3 m 


heating and cooling mode 

type 1 boundary condition. The heating region is the upper half of the 
pipe, and the cooling region is the lower half. No parasitic heat loss 
or gain is considered. 


Q = 80 W 
T = 250 K 

V 


total heat transport 
vapor temperature 



properties of working fluid and pipe material 


fluid (ammonia^ density 
fluid thermal conductivity 
aluminum thermal 'conductivity 

heat pipe operating conditions 

one-g condition 

elevation 

excess mass charge 

heat pipe performance limits 

maximum heat transport (in zero-g) 
static height 


= 670 kg/m 
= .592 W/mK 
= 180 W.mK 


ELV = .6E-2 m 
XMASS = 1.2 g. 


Q 

X 


limit 

ST 


= 215 W 
= .14E-1 m 


The sample input is gJven in Table 2. The program output for the 
above case is presented in Appendix C. The results show that the ^ 

evaporator and condenser film coefficients are .67 and 1.03 W/cm K, 
respectively, and no partial dry-out is expected. 55 iterations are 
needed to obtain the temperature distribution within .01 K accuracy. 
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NOMENCLATURE 


C's 


D 

ELV 

g 

h 


heq 

h 

par 

H 

k 

La 


Lc 

Le 

Lp 


L 

L 

N 


eff 

/ 

eff 


» 


Q 


q 


Q 


limi.'r. 


R. 

1 

R ^ 
out 


finite difference coefficients 

groove depth 

elecation 

gravitational acceleration 
over-all film coefficient 
equivalent heat transfer coefficient 

heat transfer coefficient for parasitic heat loss or gain 
puddle depth 

puddle depth at condenser end 
thermal conductivity 
length of adiabatic section 
length of condenser section 
length of evaporator section 
puddle length 

effective transport length (=La + 1/2 (Le + Lc)) 

defined in Eq. 9 

number of grooves 

total heat transport 

heat flux per unit area 

maximum heat transport 

single groove heat transport limit 

pipe inner radius 

pipe outer radius 



R 

S 

T 



w 


'^amb 


W . 
mxn 


^ST 

XMASS 


(r,ip,z) 


a 

6 

A 

y' s 
0 
P 
a 

Subscripts 

c 

e 


£ 


vapor core radius 
average lard width 
temperature 

inner surface temperature 
vapor temperature 

average outside surface temperature 
ambient temperature 
groove minimum width 
static height 
excess mass charge 

cylindrical coordinate system defined in Fig. 1 

j- 

til<" angle 

spacing between two adjacent nodes 
nodal width 

angles to specify heating and cooling zones 

puddle angle 

density 

surface tension 

condenser 

evaporator 

liquid 


w 


'wall 





TABLE 1 INPUT CARDS DESCRIPTION 


Card No. 

Format 

Name 

Description 

Unit 

Remarks 

1 

3F10.4 

XLE 

evaporator length 

m 




XLAD 

adiabatic length 

m 




XLC 

condenser length 

m 


2 

15 

NGRV 

no. of grooves 

— 



5E10.4 

GDEPTH 

groove depth 

m 




WIDTH 

average land width 

m 




ROUT 

pipe outer radius 

m 




RIN 

pipe inner radius 

m 




RV 

vapor core radius 

m 


3 

5A2 

IFLUID 

name of working fluid 

— 



10A2 

JPIPE 

name of pipe material 

— 


4 

3F10.5 

Q 

total heat transport 

W 

^ f 



TV 

vapor temperature 

K 

X J. 

KTYPE = 1 



QMAX 

max. heat transport 

W 



3F10.5 

Q 

total heat transport 

W 

f 



TCOND 

condenser surface temp. 

K 

XI. 

KTYPE = 2 



QMAX 

max. heat transport 

W 



3F10.5 

TEVP 

evaporator surface temp. 

K 




TCOND 

condenser surface temp. 

K 

if 



QMAX 

max. heat transport 

W 

KTYPE = 3 

5 

3E15.4 

CONDL 

liquid thermal conductivity 

W/mK 




RHO 

liquid density 

Kg/m^ 




CONDW 

pipe thermal conductivity 

W/mK 


6 

3E15.4 

XSM 

excess mass charge 

Kg 

Skip this card 



EL 

elevation 

~nr 

if MGRVT = 1 



XST 

static height 

m 


7 

4F10.3 

PSIE(l) 

angles to specify 

DEG. 




PSIE(2) 

heating zone 





PSIE(3) 






PSIE(4) 






Table 1 Input Cards Description (cont'd) 


Card No. Format 

8 F10.3 


9 F10.3 

E10.4 


Name 

Description 

Unit 

PSIC(l) 

PSIC(2) 

PSIC(3) 

PSIC(4) 

angles to specify 
cooling zone 

DEG. 

TEMPS 

ambient temp. 

K 2 

HPAR 

heat transfer coeff. for 
parasitic heat loss or gain 

W/m K 


Remarks 


(Repeat No. 4 to 6 cards for additional cases) 




Table 2 Sample Input 


COLUMN NO. 


000000000111111111122222222223333333333A444444AA455555555556666666666 

123456789012345678901234567890123456789012345678901234567890123456789 


INPUT DATA 


.3048 

.4572 

.1524 


27 .108E-2 .37E- 

AMMONIA ALUMINUM 

-3 .795E-2 

.544E-2 

80.0 

250.0 

215.0 


.592E 0 

.67E 3 

.18E 3 

1.2E-3^ 

.6E-2 

1.4E-2 

90.0 

270.0 

500.0 

500.0 

0.0 

0.0 

90.0 
.OE 0 

270.0 

360.0 
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EVAPORATOR 

SECTION 


ADIABATIC 

SECTION 


CONDENSER 

SECTION 



Figure 1. Heat Pipe Cross-Section 










Figure 3. Heat Pipe Model 
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Shape 3 


Shape A 



Figure 5. Puddle Shape Claasif ication 




Example 1 



= 0 (deg. ) , Ji = 60 
Ifs = 300 JV = 360 



Example 4 



ITi = 80 ,J1= 150 

= 230 , 320 


Figure 6. Heating and Cooling Zone Specification 
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PROGRAM FLOWCHART 
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V. 
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MAIN PROGRAM 
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SUBROUTINE PUDDLE 


START 


CLASSIFY 
PUDDLE SHAPE 

' 

f 

DETERMINE PUDDLE 
COVERED REGION 

' 

f 

COMPUTE TRANSPORT 
LIMIT OF EACH 
GROOVE 





RETURN 
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flS/?(SO hT Ml ft AN H 


COMPILfR OPTIONS -NAMP*' HA |N,OPT=0?fL INFCNT «B? ,S WP«nonOK , 

: soiimcf . PHr flic .NOL i st .NTioprK .lq ao, hap .mok) n , in, xkfp 

ISN nno? niMPNSioN r ( 7 1 .psh ^si .npr.psM ■^s) ,z ( 3 ? ) ,or ( 7 1 ,ops i ( 3s ) lO? ( 3?) 

li f:;. -O.Qt)3 P|HPNSJON,,JJ5t?Vj?7I,T ] ? 9 , ? 7 1 , T 1 N( ? r .2 7 i . tout I ? si. 27 I 

ISN OOOA OIMPNSION PsTf ( A I ,PSIci Al .NOPt A) .NOC ( 4 I 

JL5.N_PQD5 .UlMPNi I QN L 1 AQnO_.M ) . PS ( AODO ) ,r. T OMA ( 7 ) , mo HA ( 7 ) 

ISN OOOB niHFNSION CONOL ISI.CONOWI SI.RHOI 51 tOI SI ,TV( SI tTCOMOIS ) iTPVAPI SI , 

1 X SHI SI .H < SI . XSTI SI 

ISN nmi7 OIMPNSION UMAX ( SI .OMAXC.R ( 351 .OSINGL I 3SI 

iSN POOR OIMPNSION IPUH niS I.JPI PPI in i 

ISN OOOR COMMON C i 7 . XL P . XL A 0 1 XL C , XI. , R V , NOR V , X S T HT .L t-M , MM , N , NN , NNN , 

I MHOL . XSMA s. PL V.HMXO .OMAXr.R 

C 



C CIMP.NI = FINITP OIPPPMPNT CO F PF IT I FN1 S IW/KI 

C CONI .CONOL >= FUlin THFRMAL CO NOllC T I V IT Y (W/(M*K|| 

C f.ONW iCONQW = PI PP WALL THPRMAL CONOUCTIVITY (W/(H*KII 

C OPSTlJI « INCPPMPNT IN MSI ABOUT J (RAO. I 

C OR ( I I « INCMFMFNT in R ABOUT I (Ml 

"C 02IKI • INCRFMFMT InT ABOUT K (Ml 

r. OTSPPC » SPPCIFIPO ACCFPTABLF CHANGP IN TPMP. HPTwePN SllCCPSSIVF 

C ITPKATIOMS (Kl 

C FL.FLy » PVAPOHATOR FNO apVATION (Ml 

C FILMFV « FVAPOR'aTOR p"|L M COFFFICIFNT ( W / ( H» •2’i‘K I I 

C FILmCO - CONOPNSFR FILM CO FFFICIFNT (H/(M»»?»K|| 
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1 20X,'FLHin THERMAL CONOI IC T 1 V I T Y ! W / I M*K )) ' , T 70 , F 1 S . 5/ 

I POX. 'WAIL THERMAL CO NOI )C T 1 Vl T Y ( W / ! M»K I ) ' . T 70 . F 1 S . S ) 



ORIGIN, qr 

Of POOR Qy,,y - 


ISN 

OSS9 


WK1TF(8,913I 



ISN 

09^0 

913 

rnPMAI (////, ’9X. ‘SllRFACP ifmPFHAIIIPF OISIP IPMTIfiN 

(OFORFFS 

K 1 ' // 1 

ISN 

OSM 


UK ITF 16,9191 



ISN 

nsf.? 

919 

FOKMAl (?X, 'AX lAL ' .?0X, *f. IKCIIMFFKFNT lAL LDCAl ION' 1 



1 SN 

OSf>T 


WRITF (8,9?01 



ISN 

OShA 

930 

FFIkPAI ( IX , 'inr AT ION* 1 



ISN 

0S6S 


WRITF (8,931 1 



ISN 


931 

FnRMAT(3X,«7 ( Ml ' , ?SX, • PS 1 (OFORFFSl'/l 



ISN 

nS67 


on 9S1 J«1,HHM 



ISN 

OSAP 


nF(,PSllJ) = PSI (JI*1H0.0/PAI 



ISN 

0S69 

951 

CnNTINIIF 



ISN 

ns70 


WR I c(8,91SI (OFGPSI (JI,J = 2,MM,3) 



ISN 

0S71 

915 

FORMAT! 10X,12F9.11 



ISN 

ns7? 


WR I T F (8,932 1 



ISN 

0S7S 

932 

FORMAT! 1 



ISN 

OS 79 


no 9S2 K=?,NN,2 



ISN 

0S7S 


WRITF (8,9161 2 (K1 , ( TOUT! J,K1 , J=2,HM,31 



ISN 

ns7h 

918 

FORMAT (F7. 3, 3X , 12F9.21 



1 SN 

f)S77 

952 

CONT INIIF 



ISN 

0S7P 


WRITF (6.9B7 1 



ISN 

0S79 

987 

FORMAT (' 1 ', 35X , 'HFAT TRANSPORT OF FACH GROOVF'/I 



ISN 

nsHO 


WRITF (6.9881 (J,J=1.121 



ISN 

nsPi 


WRITF(6,9891 (OSINOL ( J1,J=2,131 



ISN 

OSH? 


WRITF (6,9881 (J,J=l?.291 



ISN 

nsfl^ 


WRITF (8, 989 1 (OS I MOL ( J1,J=19,251 



ISN 

0SP9 


WRITF (8,9881 (J.J=?S.M1 



ISN 

OSftS 


WRITF(6,9891 (OSIMOL ( JI,J=26,MM1 



ISN 

0SP6 

988 

FORMAT! // , 2X, 'OROOVF Nn.»,6X,12IBl 



ISN 

0Sfl7 

989 

FORMAT!/, 2X, 'HFAT TRANSPORT (W1',12FB.31 



ISN 

0SR8 


WRITF(6,9171 



ISN 

0S89 

917 

FORMAT! ///, 30X, 'HFAT PIPE PERFORMANCE CHA RAC TFR I S T IC S ' , // 1 


ISN 

0S90 


nTFVAP I TWFuV - TVAP 



ISN 

0S91 


OTCONn * TVAP - TWCAV 



ISN 

0S9? 


WRITF (6,9181 OTRNSP.OMAXHP, TVAP, TWFAV, TWCAV, OTFVAP 

, OTCONn, 




1 

FIL MFV, FIL MCO 



ISN 

0S93 

918 

FORMAT ( 20X, ' TOTAL HFAT TR ANSPO H T ( W 1 ' , T 70 , F 1 S . 7/ 






20X, 'MAXIMUM HFAT TR A NS PO R T ( W 1 • , T 70 , F 1 S . 2/ 






20X, 'VAPOR T FMPFKAlllRF! K 1 ' , T 70 , F 1 5 . 2 / / 






20X,'AVFRAOF FVAP. SHRFACF TFMP. (K 1 ' , T70 ,F1 5 

777 





?OX,'AVFRAOF CONO. SdRFACF T FMP . ( K 1 ' , T 70 , F 1 5 

. 2/1 





20X.'AVFRAOF FVAP. TFMPFHATIIRF PRO P (K 1 • , T70 , 
20X,'AVFRAOF CONO. TFMPFRATlIRF ORO P ( K 1 ' , T 70 , 

FIS. 2/ 
FIS.?// 





20X ,' EVAPORATOR FILM f.OF FF I C I FNT ( W / ( M*w? >»K 1 ' 

,T70,F1S. 

9/ 




?0X , 'CnNOFNSFR film f.O F F F I C I NT ( W / ( M • »K I • , T 70 . F I S . 9 I 


ISN 

0S99 


if(moravt .f(j. 11 on to 995 



ISN 

0S9(S 


!F(lPlinL .FO. 11 00 TO 966 



ISN 

nS9R 


WR ITF (6,9971 



ISN 

0S99 

99 7 

FORMAT ( /,20X, 'PlinOL F EFFECT IS INCLIIOFn'l 



ISN 

OSOfI 


00 TO 99S 



ISN 

(u,n 1 

9A8 

WR ITF (8.99H 1 



ISN 

nso? 

998 

FORMAT (/, POX , 'PlinOL F EFFECT IS NOT INCLIinFO'! 



ISN 

OhOP 

96 5 

CONT INIIF 



ISN 

0809 


no 9S3 J=2,MM 



ISN 

080S 


IFIOSINOLUI .OT. OMAXORIJII 00 TO 9S6 



ISN 

0807 

953 

CONT INIIF 



ISM 

080H 


UR 1 TF (6 ,9SS 1 




ISN nhoq P55 FflKMAI ( /,?OX , 'Ml PftkIIAL OKY-mn IS FXPFCTFn't 

IbN Ohio r.n Tfl 9SA 



e?SiiCfiV... t-i 

OF POOK QUAuskV 


ISn OMl 

LiN._Qiiii’„. 

ISN n^.n 

ISM nM4 

NSS WM n 1 ^ ( 6 ,qs7 1 

NS7 PCIKMAI ( / , ^()X, 'PAH 1 1 m OHY-llUJ J 

9S6 ( ONT INIIP 

W« n P ( (S . 97 1 1 Ml IN 

S_.fr XPFf 1 FOj 1 


ISN (IMS 
ISN QA16 , 

971 FnwMAI ( //,20X. ' iniAL NllMWeK flF 
_ HQ TH SO . 

ITFRAT IONS RF(.i|l 1 Rho * ' 1 

IS) 

ISM QM 7 

c 

97^ . f.nisiT INllF.. . . 



ISN niSl fl 
ISM fiAi q 

WRITF(6.961I IHIIN 
WH 1 IF (h-97S 1 



iSN 06?0 

976 FOKMA 1 ( // , 20X , • ini AL HFAT IMPlII 

FXCFFOS MEAT TRANSK1HT 

LIMIT ' 1 

ISN, 0671 . 

C 

SO continue 

. 


ISN fl6?? 
ISN Qh?-^ 

STOP 

FNI1 





ORSGfiWiL ?i:cs 
OF POOR Ql-r^liJY 


4^ 


(■(‘.'/•■‘fin HlWIKAN H 


COMP ILFR hpT iriMS -”namf= ma iNinPiVo?,! Tmfcnt = h?~.s r?F = finnnK , ' 
, sniiPCF. Furmr. .Nf^L 1ST .I'ftoFCK .Lnfto.MftP.^nhpn , iii.xhff 


ISN 

nno? 


SlihRmniNF PlinOLF 



c 

L. 

THIS SllHROiniNF OFTFRMINFS THF RFOION f.nwFHFO HY A 1, Inillll PIKIL ANil 

iMpnsFS cnNoniOM that hfat ipamsffr ratf is ;fpo hvfm ihf RFf.iiiN 



c 

r 

THIS PROr.HAM ALSO CALULATFS CHAMf.F OF MAX. HFA1 TRANSPORT OOF TO 
pimni F FOR faCh ohoopf 

I SM 

nonx 

c 

n IMFNSinN Cl AOOO ,R l . ? I l .OMAXORI PSI , XPOL I PF I , XTWNS ( i 

ISN 

nooA 


COMMON C . 7 , XL F . XL AO.XLf f XL , R V , NOR V , X S T H T .L , M , MM , N , NN , NNN . 
1 RHOL .XSMAS, Ft V , 1 tMXO , OP A XOR 

f. 

c 

ISN 

ISN 

c c 
c c 
c c 


MPARII.J.Kl = (K-n*L*M ♦ (I-n*M ♦ J 
PAI = 2.1AIS4 

ISN 

0007 

c 

on 740 J»?,MM 

I SN 
1 Sm 

000 P 
OOOP 


XPOI ( J 1 ■ 0.0 
740 (,OnTIN(/F 

ISN 

0010 

f, 

XSVOL « XSMAS/RHOl 



c 

CLASSIFY PtinOLF SHAPF 


Tsn noTl ’ '"xHr~= 2 . 0 *)(SVnL /( PaI »R V»XL) 'i , 5 *FLV 

I^M no I? KHO = XHl - FLV 


ISM oon iFiFuv .r,i. ( 2 .o*rvm cn to 702 

IS M OQIS I -F-IXM O. .i f. O.P) r. n lo 710 


ISN 

0017 

IF (XHl 

.OT. 

(?.0*RV1I GO TO 

71 1 





1 SN 

noi<5 

on TO 712 







ISN 

0020 

703 CONTINIIF 







1 SN 

0021 

IF IXHl 

• 1 F. 

I2.0»RVM on TO 

710 





ISN 

002 T 

IF ( XHO 

.07. 

0.01 on TO 711 






I SN 

002 S 

on TO 713 








c 

ISM flOPh 710 CONT IMIIP 


c, 

I PliPniE -5H.A.Ff 1 


ISN 

0027 

C 

XHl ■= 2 . 0 *RV>S 0 R T ( XSVOL *a V/ ( PA I*XL» I P VI ••31 1 


ISN 
I Sn 

00 2 R 
00 2 R 


XLP « XHMXL/FLV 
I ■= ? 


ISN 
ISN.. 

0030 

0031 


on 770 KK=l ,N 

K = ^^M^l - KK 


ISN 
I SN 

00 3 2 
00^3 


on 720 J= 2 .MM 

MP = MP AH ( 1 , J, K ) 


ISN 

ISN 

0034 

003 5 


7 HC = XL - Z( X 1 

lF( 7 Hf. .0 1 . XI PI GO TO 721 


ISN 

JS'-I 

00 ■^7 
OO'^P 


HG,P = 1.0 - XHl/RV + 7 Hr.*FL V/(RV*XL ) 
IF(HG,P ,G, 1 . 1 . 0 ) Hr.P= 1.0 


ISN 
1 SN 

0040 
004 1 


PANOL = ARCnSIHOPI 

.IPHO = P ANG.L »P r.R V/ ( 2 . “PA 1 ) 


I Sn 

LSN 

OOA? 

oriA 4 


IFU . 01 . UPHO +11 ,aho. j .LI. (MM-jpnn+in fo to 72 o 
r (Mp,j ) = 0.0 


ISM 

1 SM 

004 F 
(1040 

7 20 

XPIH ( J 1 * ?Hf. 

rnta inhf 


ISN 

ISN 

004 7 
004 H 

721 

CONT INIIF 

r.n TO 7 F 0 
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I SN 

0049 

■ r. 

71 1 

f.ONT INIIF 




c 

r. 

RiinoLF Sharp ? 


I SN 

OOSO 

c 


XLPl = XL - SnPTt4.0*XL*(RAI*XI.*(RVI«4? - XSVXIL 

1/ ( rA I *hu*FLV I 1 

ISN 
1 SN 

OOSl 

OOS? 



XHO = 2.0*RV - IXL - XLR1I»FLV/XL 
I = ? 


ISN 
I SN 

OOS^ 

0044 



on 7?S KKS 1 ,N 

K = MNN - KK 


ISN 
I SN 

OOSS 

flOSiS 



on 7?S J=?*MM 

MP c HRAP ( I , J, K 1 


ISN 
1 SN 

00S7 

OOSfl 



7HC = XL - 2(Kt 

lP(?Hf. .01. xLRii on in 7?o 


ISN 
1 SN 

0000 

OOM 



r.iMP.n = 0.0 

OONT INHF 


ISN 
I SN 

000? 

000^ 



PANOL = APCnS((7Hf. - XL PI l*FL U/( RV»XL 1 - 1.01 
JROO = PANOL *NOR V/ 1 ?. OPA I 1 


ISN 

ISN 

0004 

0000 



1F(J .01. (aRit04ii .ano. j .LI. I MM-.1PU0+1 1 1 rn 
C (MP. 1 1 « 0.0 

in 72 s 

ISN 
I SN 

0007 

OOOP 


7?S 

XPOLUl = 2HC 
CONT INIIF 




ORIGffMAL PAGS jg 
OF POOR QUALITY 


1 SN 

nnsq 

on TO 7S0 
c 



ISN 

n070 

71 ? CONT INI IF 
f 





C , PIIOOLE SHAPF 3 

r. 



1 SN 
ISN 

no7i 

no7? 

XHl « ?.0*XSva /(PAI 1 ♦ .S*FLV 

I • ? 



ISN 

iSN 

nri7-» 

007* 

on 7S0 KKsl.N 
K = NNN - KK 



ISN 007S 
LSN 007f. 

on 730 J=?,MM 
MP = NPAR (1 1 J . X * 



ISN 

ISN 

0077 

007fi 

7HC « XL - ?(KI 

PANOI . APfnSIl.O - XHl/HV ■» ?HC»FL V/ (RV*XI. II 



ISN 

ISN 

0079 

OOHO 

JPUn « P ANOL *MOR V/ ( ? . *PA I 1 

IFU .01. tJPimrll .ano. J .LT. (mh-jpiio+ 1 I I m 

TO 

730 

ISN 

ISN 

OOIi? 

OOH'^ 

f, IMP, 1 1 « 0.0 
XPIILIJI ' 2HC 



ISN 

ISN 

OOR* 

OOHS 

730 CONT INIIF 
on in 7S0 



HPI 

ipm 

r. 

713 CONT INIIF 





C 

r. PIIOOLE SHAPF * 




mtitt 

C 

XlPl « XL - SOP 1 ( *. 0»XI • ( PAI»X1, • ( «VI ••? - XSWOL 1 / ( PA| •M V*FL VI 1 

ISN 

ISN 

OORR 

00H9 

XLP » XLPl ♦ ? .0»P V*XL /EL V 
1 • ? 



I SN 
ISN 

0090 
009 I 

on 73S KK*1-N 
K = NNN - KK 



1 SN 
ISM 

009? 

009-^ 

on 73S J=?,MM 
MP = HPAR ( I , J . K 1 



I SN 
ISN 

009* 

0099 

7HC = XL - ZIKI 

IF ( 7HC .01 . XI P 1 1 on in 736 



I SN 
ISN 

0097 

009R 

f. (MP. 1 1 « 0.0 
73 6 rflNTINIlF 



I SN 
ISN 

0099 

0101 

iF(7Hf, .0 1. XL PI on in 737 

PANOL = ARCnSKZHC - XL PI 1 *FL V/ ( H V* XL 1 - 1.01 



ISN 

ISN 

010? 

OlOR 

JPIIO = PANOL •NORV/(?.*PA| 1 

iFi.i .01. (jPiin+11 .ANO. J .LT. 1 MM-jpiin+1 1 1 m 

1(1 

735 

1 SN 
ISN 

OIOS 

OlOiS 

C. (MP. 1 1 « 0.0 
XPDUJI = ZHC. 



I SN 
ISN 

0 107 
0 1 OR 

73S CONT INIIF 
73 7 fflNl INIIF 



ISN 

0109 

C 

7S0 CONT INIIF 



ISN 

Olio 

C 

on 7*1 J=?.MM 



1 SN 
ISN 

0111 

Oils 

IFIXPOLUI .01. XL Cl on TO 7*? 

XTRMSUl = .S*XLF + XLAn ♦ .S*(XLC - XPOLUIl 



I SN 
ISN 

Oil* 

Oils 

on Tn 7*3 

7*? IFIXPOLUI .01. (XLC+XLAntI On Tn 7** 



I SN 
ISN 

0117 

OUR 

XTKNSUl « .S*XLF + (XLC XLAO - XPOLUIl 

on Tn 7*3 



I SN 
ISM 

0119 
0 1 ?0 

7** XTHNSUI - .S*(XL - XPOLUIl 
7*3 CONT INIIF 



ISN 

ISN 

01 ?1 
01 ?? 

XTRNSO = .5*XLF ♦ XLAO ♦ .5*XLC 

OMAXORUl e OMXO* (XTRNSO/XIRMSUI 1* 1 1 .0 - IFLV/XSIHII* 
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SAMPLE OUTPUT 



V 












THFkHAl. ANALYSIS 

OF AXIALLY GROOVED HFAT PIPF 
















MAMFS OF WORKING FLUIO AMO PIPF 

MATFRIAL 







working fliiio 
PIPF MATFRIAL 

AMMONIA 
AL IIM I NIIM 






HFAT PIPF OIMFNSIflN.S 



FVAPORATOR LFNGIHIMI 

U.^IIAHHF 01) 







AOIAHATir. LFNGTHIMI 
f.OMOFNSFR LFNGIHIMI 

0.4S7P0F 01) 
0.1S2A0F 00 







TOTAL PIPF LFNGTHIMI 
PIPF OHTFH OIAmFIFRImi 

0.91AA0F 00 
0. 1 S9H0F-O1 







PIPF INNFR OIAMFIFRIMI 
VAPOR r.OHF HAOlllSIMl 

0. lOHhOF-01 
0. ASOODF-02 






fiHfinVF niMENSIONS 



NIIMHFH OF GROOVES 
GROOVE OFPTHIMI 

?7 

0. lOHIU)F-()2 







AVFRAGF LANO WIOTHIHI 

0.37000F-03 





MFATING AMO CnriLlMG hoofs 



EVAPORATOR RFGION 

ATM-UNIFORM HEATING 








HEATING RFGION COVFRS FROM PSI= 90.0 

OEG 

TO 

PS 1 = 270.0 

OEG 



CONOENSFR RFGION 

NON-IINIFORM COOLING 

COOLING RFGION COVFRS FROM PSI= 0.0 

OFG 

Til 

PSI= 90.0 

OFG 




Amo FROM PS1=270.0 

DEft 

T7T 

■psi='36ft'.n 

"UFn 


o o 

™ S3 



O 

33 n 

3> '-i 


r • 






Cask NiiMii(-'K 


1 


( 




HhiU.ai.li’.f l l f'Hlimhli Lllf' t ll T lliN^a 


ONF-c. f nNO n iiiN 
J YPlEi... t linuilDAKy-i.LlA'DLlI[Ifl 


TriT Al HF A 1 I Ml^lll ( W I hll . (1(1 


VAPHR TFMPFRATIIKFIKl 
AMhlFMT Itl-'HFKATliHFtK ) 

2S(l.(lO 
(1 .0 




mFAT TMANSFFH CflFFF. FUR PARASITIC 
HFA? LOSS OH r.AIOlW/l M**?*K) 1 

o.u 

— - — - - 



Fxr.FSS MASS CHAHCFIKO) 
FLFVAT mO(M) 

U. 120flliF-fi2 
(1. h()nnnF-ii2 



STATIC HFir.HT OF PIPFIMl 

o.uonoF-ni 




PROPFHTIFS OF WORKING FL Ul 0 AND PIPF HATFRIAL 

FLUID riFMSI 1 Y (KG/M»»^1 

0.670(IUF (IT 




FLiiin thermal cnNOoCTi viTy(w/iM*Ki I 

O.S420nF 00 




WALl THERMAL CO NOl iC T I U I T Y I H / ( M» K 11 

0. IHOnOF 03 




• 


SllRFACF UMPERATURF 0 1 SI R 1 BUT in N (nFGKtFS K) 


AXIAL 



C IKCtlMFFRFNTl AL LOCATION 





Lnr.AT ION 
7 (HI 




PSI 

intGKFFSl 







6.7 

56.7 

«6.7 

126.7 

lftft.7 

206.7 

256.7 

286,7 

326.7 


O.OIS 

250.77 

250.95 

251.55 

251.80 

251.93 

251 .87 

?51 .60 

251 .07 

250.81 


0.07A 

0.137 

250.78 

250.78 

250.96 

2<^0.'7h 

251.55 
. 251.55 

251.81 
251 .81 

251.95 
251 .95 

251.88 
251 .88 

251 .62 
251.62 

251 .08 
251.08 

250.82 

250.82 

- 

n.iQfl 
n. 259 

250.78 
250. 7H 

250.9ft 

250.95 

251.55 
251 . 55 

251.81 
251 .81 

251.95 

251.95 

251.88 
251 .88 

251.62 
251 .61 

251.08 

251.08 

750.82 

250.82 


0.351 

0.533 

250.02 

250.00 

250.02 

250.00 

250.02 

250.00 

250.03 

250.00 

250.03 

250.00 

250.03 
2 50 .00 

250.03 

250.00 

250.02 

250.00 

250.02 

250.00 


0.71ft 

0.7R5 

259.98 

256.81 

259.99 

257.61 

2A9.99 
2AH . 7ft 

250.00 

259.52 

250.00 

259.63 

250.00 

259.55 

259.99 
259.1 1 

259.99 

258.06 

259.99 

257.22 


0.H15 

O.hAft 

2A6.1R 
i •*! 7 

257.12 

257.11 

258.56 

258.55 

259.33 
259. 33 

259.58 

259.58 

259.59 

259.59 

259.01 

259.00 

257.79 

257.78 

256.60 

256.59 


0.b7ft 

Oi'jgT 

256.17 

256.17 

257.11 

257.11 

258.55 

258.55 

2 59 .'3 3 
2 59 . 33 

259.58 

259.58 

259.59 

^59.59 

259.00 

259.00 

257.78 

257.78 

256.59 

256.59 
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HHAI lt<A^.SHl^<T fit- l-Af.H {;t<nnVF 


GRnnve nd. 


1 

2 

3 

A 

5 

6 

7 

a 

4 

“TO 

11 

— n 

HFAT TRANSPORT 

(W) 

1 ,7n 

1 .7H5 

1.911 

2.09A 

2 . 3A0 

2 .652 

3.015 

3.331 

3.5H5 

3.7H 1 

3.92A 

A.ni H 

f.RfinvF NO. 


13 

LA . 

15 

16 

17 

1 H 

1 9 

20 

21 

22 

2 3 

2a 

HFAT TRANSPORT 

(M) 

A .067 

A , 069 

A. 027 

3.9 39 

3.H00 

3.60H 

3.356 

3.0A2 

2.67R 

2 . 3 6A 

2.115 

1.927 


fifionvF NO. 


25 

26 

27 










HFAT TRANSPORT 

(H) 

1.796 

1.71U 

1.690 











HFAT P'lPF P'FRFOR'MANCF CHARAC'tEM IST ICS' 


total HFAT transport (W) 

78. 3A 


maximum HFAT TRANSPORT(W) 
VAPOR TFMPFRATURFtKI 

132.15 

250.00 


AVFRAGF FVAP. SdftFAf.F TFMP.IKI 

251.36 


AVFRAGF CONO. SliRFACF TFMP.IKI 

2A8.23 


AVFRAGF FVAP. TFMPFRATljRF OROP(K) 

1 .36 


AVFRAGF CONO. TFMPFRATURE OROPIKl 

1.77 


EVAPORATOR FILM C OF F F 1 f. I FN T I W / I M* » 2»K ) 

0.6701F Oa 


CONOENSFR FILM COFFF 1 C 1 NT ( W/ 1 M**2»K 1 

0.102HF 05 


PlinOLE EFFECT IS INCLUDED 

NO PARTIAL DRY-nilT IS EXPECTED 

TOTAL NIIMHFR OF ITERATIONS RFOUIRFD = 55 



o , 


<Q 

C 

> S 

-< Cft 



